

use "HS_Inequality_CS.dta", clear

*Run Soloman Hsiang's Ado file for OLS Spatial HAC SEs
run "ols_spatial_HAC.do", nostop


*Determine minimal sample containing all relevant variables to ensure stable samples for tables
quietly eststo: xi: ols_spatial_HAC gini1935 pctho  res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi , lat(ycen) lon(xcen) timevar(yar) panelvar(GISJOIN) dropvar distcutoff(250)

gen Samp1 = (e(sample)==1)


quietly  xi: ols_spatial_HAC GiniAg1935 pctho  res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi , lat(ycen) lon(xcen) timevar(yar) panelvar(GISJOIN) dropvar distcutoff(250)

gen Samp2 = (e(sample)==1)

*Appendix Table A1 (Cross Section Summary Statistics)
********************************************************************************
********************************************************************************
********************************************************************************
eststo clear
estpost sum gini1900 gini1935 GiniAg1900 GiniAg1935 RD pd_1860 pcthome elevationMean roughness precipMean tempMean PctHighsoil StreamDensity   RailDensity  res2010 if Samp2 ==1

*still refining this...
esttab , cells("count mean sd min max")  nomtitle nonumber noobs




*Table 1 and Table A3: The Relationship Between Homesteading and Land Inequality 
********************************************************************************
********************************************************************************
********************************************************************************
eststo clear

capture quietly xi: ols_spatial_HAC lngini1935 pctho  i.decade i.state constant if Samp1 ==1, lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_1

quietly xi: ols_spatial_HAC lngini1935 pctho   elevationMean roughness precipMean tempMean PctHighsoil StreamDensity   RailDensity  res2010 i.decade i.state constant  if Samp1 ==1, lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_2

quietly xi: ols_spatial_HAC lngini1935 pctho  res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if Samp1 ==1 , lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_3

quietly xi: ols_spatial_HAC lnGiniAg1935 pctho  i.decade i.state constant if Samp2 ==1, lat(ycen) lon(xcen) timevar(yar) panelvar(GISJOIN) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_4

quietly xi: ols_spatial_HAC lnGiniAg1935 pctho   elevationMean roughness precipMean tempMean PctHighsoil StreamDensity   RailDensity res2010   i.decade i.state constant  if Samp2 ==1 , lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_5

quietly xi: ols_spatial_HAC lnGiniAg1935 pctho  res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi   if Samp2 ==1 , lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1935 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto T1_OLS_6

esttab T1_OLS_*,  label  indicate("State FE = _Is*" "Arrival Decade FE = _Id*") se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps keep(pcth* elevationMean roughness precipMean tempMean PctHighsoil StreamDensity   RailDensity  res2010) scalar(df_a_initial MDV)




*Table 3: Homesteading along the 1860 Frontier
********************************************************************************
********************************************************************************
********************************************************************************
local i= 1

foreach D in  "200000" "150000" "100000" "75000" "50000" "0"{

quietly xi: ols_spatial_HAC  pctho RD  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_f <= `D',  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum pctho if e(sample) ==1
estadd scalar MDV =r(mean)
est sto HS_OLS_RD_`i'

local i = `i'+1

}
esttab HS_OLS_RD_*,  label  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps keep ( RD   ) title(Effect of 1860 Frontier on Homesteading) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar( MDV)





*Table 4: Homesteading and Historical Land Inequality along the 1860 Frontier
********************************************************************************
********************************************************************************
********************************************************************************
local i =1
*Table and Appendix Tables: Effect of Homesteading on Farm Size Gini Along the Frontier
foreach D in  "200000" "150000" "100000" "75000" "50000" "0" {

quietly xi: ols_spatial_HAC lngini1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto gRD_1900_`i'


quietly xi: ols_spatial_HAC lnGiniAg1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum GiniAg1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto G_Ag_RD_1900_`i'

local i = `i'+1
}

esttab gRD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title(Effect of  Homesteading on 1900 Land Endowment Gini) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)

esttab G_Ag_RD_1900*,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title(Effect of  Homesteading on 1900 Farm Size Gini) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)



*Table A4: Homesteading along the 1860 Frontier, Omitting Western States
********************************************************************************
********************************************************************************
********************************************************************************
local i= 1

foreach D in  "200000" "150000" "100000" "75000" "50000" "0"{

quietly xi: ols_spatial_HAC  pctho RD  res2010  i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_f <= `D' & West ==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum pctho if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nwHS_OLS_RD_`i'



local i = `i'+1

}

esttab nwHS_OLS_RD_*,  label  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps keep ( RD   ) title(Effect of 1860 Frontier on Homesteading) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)


*Table A5: Homesteading and Historical Land Inequality along the 1860 Frontier, Omitting Western States
********************************************************************************
********************************************************************************
********************************************************************************

local i =1
*Table and Appendix Tables: Effect of Homesteading on Farm Size Gini Along the Frontier
foreach D in  "200000" "150000" "100000" "75000" "50000" "0" {

*No West version

quietly xi: ols_spatial_HAC lngini1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nwgRD_1900_`i'

quietly xi: ols_spatial_HAC lnGiniAg1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum GiniAg1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nwG_Ag_RD_1900_`i'

local i = `i'+1

}

esttab nwgRD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Land Endowment Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)

esttab nwG_Ag_RD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Farm Size Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)






*Table A6: Homesteading along the 1860 Frontier, Omitting Oklahoma
********************************************************************************
********************************************************************************
********************************************************************************
local i= 1

foreach D in  "200000" "150000" "100000" "75000" "50000" "0"{

quietly xi: ols_spatial_HAC  pctho RD  res2010  i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_f <= `D' & West ==0 & state !="OK",  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum pctho if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nookHS_OLS_RD_`i'



local i = `i'+1

}

esttab nookHS_OLS_RD_*,  label  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps keep ( RD   ) title(Effect of 1860 Frontier on Homesteading) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)



*Table A7: Homesteading and Historical Land Inequality along the 1860 Frontier, Omitting Oklahoma
********************************************************************************
********************************************************************************
********************************************************************************
local i =1
*Table and Appendix Tables: Effect of Homesteading on Farm Size Gini Along the Frontier
foreach D in  "200000" "150000" "100000" "75000" "50000" "0" {

*No West version

quietly xi: ols_spatial_HAC lngini1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0 & state !="OK",  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nookgRD_1900_`i'

quietly xi: ols_spatial_HAC lnGiniAg1900 pctho  res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0 & state !="OK",  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum GiniAg1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto nookG_Ag_RD_1900_`i'

local i = `i'+1

}

esttab nookgRD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Land Endowment Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)

esttab nookG_Ag_RD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Farm Size Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)






*Table A8: Homesteading along the 1860 Frontier, Omitting Reservation Counites
********************************************************************************
********************************************************************************
********************************************************************************
local i= 1

foreach D in  "200000" "150000" "100000" "75000" "50000" "0"{

quietly xi: ols_spatial_HAC  pctho RD    i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_f <= `D' & West ==0 & res2010==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum pctho if e(sample) ==1
estadd scalar MDV =r(mean)
est sto norezHS_OLS_RD_`i'



local i = `i'+1

}

esttab norezHS_OLS_RD_*,  label  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps keep ( RD   ) title(Effect of 1860 Frontier on Homesteading) mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)




*Table A9: Homesteading and Historical Land Inequality along the 1860 Frontier, Omitting Reservation Counties
********************************************************************************
********************************************************************************
********************************************************************************
local i =1
*Table and Appendix Tables: Effect of Homesteading on Farm Size Gini Along the Frontier
foreach D in  "200000" "150000" "100000" "75000" "50000" "0" {

*No West version

quietly xi: ols_spatial_HAC lngini1900 pctho     i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0 & res2010==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum gini1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto norezgRD_1900_`i'

quietly xi: ols_spatial_HAC lnGiniAg1900 pctho     i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  if abs_f < = `D' & West ==0 & res2010==0,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
quietly sum GiniAg1900 if e(sample) ==1
estadd scalar MDV =r(mean)
est sto norezG_Ag_RD_1900_`i'

local i = `i'+1

}

esttab norezgRD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Land Endowment Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)

esttab norezG_Ag_RD_1900* ,  label keep ( pctho*  )  indicate("State FE = _Is*" ) se(a3) b(a3) star(* 0.1 ** 0.05 *** 0.01) ar2  replace  nogaps  title("Effect of  Homesteading on 1900 Farm Size Gini (No West)") mtitles( "200km" "150km" "100km" "75km" "50km" "0 km") scalar(df_a_initial MDV)








*Figure 3 & Appendix Figure A7
********************************************************************************
********************************************************************************
********************************************************************************


foreach t in  "1910" "1920" "1925" "1930"  "1940" "1950"  "1978" "1982" "1987" "1992" "1997" "2002" "2007" "2012" "2017"{

gen lnGiniAg`t' = ln(GiniAg`t')
}




local i = 1
foreach t in "1900" "1910" "1920" "1925" "1930" "1935" "1940" "1950"  "1978" "1982" "1987" "1992" "1997" "2002" "2007" "2012" "2017"{


quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi  ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_v1

quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010   i.decade i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if Samp2 ==1 , lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_v2

quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010    i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=200000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_200k


quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=150000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_150k

quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010   i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=100000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_100k

quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010    i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=50000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_50k


quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010    i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=75000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_75k


quietly xi: ols_spatial_HAC lnGiniAg`t' pctho res2010    i.state constant i.pctRugg i.pctPre i.pctElev i.pctTem i.pctStre i.pctRa i.pctSoi if abs_ <=000 ,  lat(ycen) lon(xcen) timevar(yar) panelvar(fips) dropvar distcutoff(150)
est sto olsgini_plot`i'_0k

local i = `i'+ 1

}







*Figure 3
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_v2, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_v2, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_v2, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_v2, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_v2, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_v2, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_v2, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_v2, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_v2, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_v2, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_v2, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_v2, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_v2, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_v2, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_v2, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_v2, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_v2, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr




*Figure A7, 200 km Sample
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_200k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_200k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_200k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_200k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_200k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_200k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_200k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_200k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_200k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_200k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_200k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_200k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_200k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_200k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_200k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_200k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_200k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr


*Figure A7, 150 km Sample
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_150k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_150k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_150k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_150k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_150k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_150k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_150k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_150k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_150k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_150k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_150k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_150k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_150k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_150k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_150k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_150k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_150k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr

*Figure A7, 100 km Sample
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_100k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_100k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_100k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_100k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_100k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_100k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_100k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_100k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_100k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_100k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_100k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_100k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_100k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_100k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_100k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_100k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_100k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr


*Figure A7, 75 km Sample
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_75k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_75k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_75k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_75k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_75k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_75k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_75k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_75k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_75k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_75k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_75k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_75k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_75k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_75k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_75k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_75k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_75k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr


*Figure A7, 50 km Sample
********************************************************************************
#delimit ; 
coefplot (olsgini_plot1_50k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_50k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_50k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_50k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_50k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_50k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_50k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_50k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_50k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_50k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_50k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_50k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_50k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_50k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_50k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_50k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_50k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr



*Figure A7, Border Sample
********************************************************************************

#delimit ; 
coefplot (olsgini_plot1_0k, offset(0.1) mcolor(gs5) ciopt(color(gs5)) msymbol(circle))
(olsgini_plot2_0k, offset(02) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot3_0k, offset(04) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot4_0k, offset(05) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot5_0k, offset(06) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot6_0k, offset(07) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot7_0k, offset(08) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot8_0k, offset(010) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot9_0k, offset(012.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot10_0k, offset(013.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot11_0k, offset(014.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot12_0k, offset(015.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot13_0k, offset(016.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot14_0k, offset(017.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot15_0k, offset(018.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle))
(olsgini_plot16_0k, offset(019.5) mcolor(gs5) ciopt(color(gs5))  msymbol(circle) )
(olsgini_plot17_0k, offset(020.5) mcolor(gs5) ciopt(color(gs5)) msymbol(circle) )
, vertical keep(pctho*) yline(0, lpattern(dash) lcolor(red)) legend(off) scheme(s1mono)   ytitle("{&part}Gini (Farm Size)/{&part} % Homesteads") 
xlabel(0 " " 1 "1900" 3"1910"  5 "1920" 7 "1930" 9 "1940" 11 "1950"  13.5 "1978" 15.5 "1987" 17.5 "1997" 19.5 "2007" 21.5  "2017" 22.5 " ") ;
#delimit cr

